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PAIRING OF SUPERMASSIVE BLACK HOLES IN UNEQUAL-MASS GALAXY MERGERS 
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ABSTRACT 

We examine the pairing process of supermassive black holes (SMBHs) down to scales of 20-100 pc using 
a set of Af-body/SPH simulations of binary mergers of disk galaxies with mass ratios of 1:4 and 1:10. Our 
numerical experiments are designed to represent merger events occurring at various cosmic epochs. The initial 
conditions of the encounters are consistent with the ACDM paradigm of structure formation, and the simula- 
tions include the effects of radiative cooling, star formation, and supernovae feedback. We find that the pairing 
of SMBHs depends sensitively on the amount of baryonic mass preserved in the center of the companion 
galaxies during the last phases of the merger. In particular, due to the combination of gasdynamics and star 
formation, we find that a pair of SMBHs can form efficiently in 1:10 minor mergers, provided that galaxies are 
relatively gas-rich (gas fractions of 30% of the disk mass) and that the mergers occur at relatively high redshift 
(z ~ 3), when dynamical friction timescales are shorter. Since 1:10 mergers are most common events during the 
assembly of galaxies, and mergers are more frequent at high redshift when galaxies are also more gas-rich, our 
results have positive implications for future gravitational wave experiments such as the Laser Interferometer 
Space Antenna. 

Subject headings: black hole physics — cosmology: theory — galaxies: interactions — hydrodynamics — 
methods: numerical 



1. INTRODUCTION 

Compelling dynamical evidence indicates that supermas- 
sive black holes (SMBHs) with masses ranging from 10 6 to 
abov e 10 9 M w reside at the ce nters of most galactic spheroids 
(e.g.. iFerrarese & Ford! 12005). The masses of SMBHs cor- 
relate with various properties of their hosts, e.g. lum inos- 
ity or mass ( Magorria n& al.lll998t iHaring &Rixll2004l) and 
veloc ity dispersion (Ferrarese & Merritt 2000; Gebhar dt & alj 
2000). In the currently favored model for structure formation, 
the ACDM cosmology, galaxies grow hierarc hically through 
merg ers and accretion of smaller systems (e.g JWhite & Reei 
1978). Thus, if more than one of the merging galaxies con- 
tained a SMBH, the presence of two or more SMBHs in their 
merger remnant will be almos t inevitable during galaxy as- 
sembly (Begel man et al.| [l98Q). However, it is unclear if the 
dynamical processes at play are efficient in forming a close 
SMBH pair with separations ~ 10-100 pc, which may subse- 
quently shrink to a bound binary, and eventually merge via 
gravitational wave radiation. Such black hole coalescence 
events are expected to give rise to gravitational wave bursts 
that should be d etectable by th e Laser Interferometer Space 
Antenna (LISA) dVecchioll2004l) . 

SMBH pairing has been shown to proceed quickly 
when both compact objects are hosted by steep stel- 



Electronic address: callegar@physik.uzh.ch 

1 Institute for Theoretical Physics, University of Zurich, Winterthur- 
erstrasse 190, CH-9057 Zurich, Switzerland 

■ Institut fiir Astronomie, ETH Zurich-Honggerberg, Wolfgang-Pauli- 
Strasse 16, CH-8093 Zurich, Switzerland 

3 Center for Cosmology and AstroParticle Physics; and Department of 
Physics; and Department of Astronomy, The Ohio State University, 191 West 
Woodruff Avenue, Columbus, OH 432 1 US A 

4 Universita Milano Bicocca, Dipartimento di Fisica G. Occhialini, Piazza 
della Scienza 3, 1-20126 Milano, Italy 

5 Department of Astronomy, University of Washington, Box 351580, Seat- 
tle, WA 98195, USA 

6 Department of Physics and Astronomy, McMaster University, 1280 Main 
Street West, Hamilton, ON L8S 4M 1 , Canada 



lar cusps approaching each other from close distances 
(Milosavli evic & Merrittll200l1) . or when embedded in a cir- 
cumnuclear gaseous disk unde r appropriate thermodynamic 
conditions ( Ma yer et al J 12007). but whether the large-scale 
merger can lead the SMBHs to such a favorable configuration 
is still a matter of debate. Previous studies found that, fol- 
lowing a galaxy merger, the relative distance of the SMBHs 
in the remnant is very sensitive to th e structure of the merg - 
ing galaxies, and t o their initial orbit dGovernato et al.lfl994h . 
Kazantzidis & al. (2005) showed that pairing is efficient in 
equal-mass disk galaxy mergers with cosmologically relevant 
orbits, while the presence of a dissipative component is nec- 
essary for the pairing of SMBHs in 1 :4 mergers. Other re cent 
studies (e.g.. ISpringel et all 120051 Uohansson et al.ll2009h fo- 
cused on the effect of energetic feedback from black hole ac- 
cretion onto the surrounding galaxy, but were not designed to 
follow the orbital evolution of SMBHs. Substantially less ef- 
fort has been devote d to examining the fate of SM BHs in mi- 
nor mergers (but see iBovlan-Kolchin & Ma 2007) , which are 
much more frequent in AC DM cosmologies (lLacey & Cole! 
ll993tlFakhouri & Ma 2008). Investigating the necessary con- 
ditions for SMBH pair formation in this regime is of pri- 
mary importance for the search of gravitational waves and for 
SMBH demographics and activity. 

In this Letter, we report on the efficiency of the SMBH pair- 
ing process using a set of A^-body/SPH simulations of disk 
galaxy mergers, with mass ratios q = 0.25 and 0.1, constructed 
to represent mergers occurring at various cosmic epochs. The 
choice of the initial conditions, in particular the masses of the 
SMBHs, is such that the corresponding SMBH coalesce nce 
events would be detectable with LISA (Ses ana et al.l l2005). 



2 



Callegari et al. 



2. SIMULATION SET-UP 

The galaxy models were initialized as thr ee-component sys - 
tems following the methodology outlined in lHernquist] (1 1993). 
They comprise a Hernquist spherical stellar bulge (Hernquist 
1990), an exponential disk with a gas mass fr action f e , and an 
adiabatically contracted dark matte r halo jB lumenth al et al.l 
119861) with an initial NFW profile dNavarro et alj|1996l) . A 
collisionless particle representing the SMBH was added at the 
center of each galaxy. 

Our reference model is a Milky-Way type galaxy, with a 
virial velocity V v i r = 145 km/s, a disk mass fraction Mj = 
0.04M vir , and a bulge mass fraction M b = 0.008M vil . The 
mass of its central SMBH is M Bil = 2.7 x 10 6 M n , consis- 
tent w ith the updated M BH -Mb u ig e relation (Harin g~& Rixl 
2004). The disk scale height and the bulge scale radius are 
Zo = 0-lRd and a = 0.2/?,/ respectively, Rj being the exponen- 
tial disk scale length . Rj is determined following the model by 
iMo. Mao & White! (1 19981) (MMW hereafter), which yields 
disk galaxies lying on the Tully-Fisher relation. Models at 
redshift z = w ere initialized with a halo concentration pa- 
rameter c = 12 (Bul lock et al.l l2Q0~Tb . We also ran mergers 
with initial conditions rescaled to z = 3 according to MMW, 
keeping V v j r fixed, as ex pected for the progenitors of our z = 
models (iLi et al.ll2007t) . Considering high-redshift mergers 
is crucial, because the merger rate increases with look-back 
time, and a large fraction of the gravitational wave signal 
from coalescences of SMBH binaries is predicted to origi- 
nate from this cosmic epoch at the corre sponding mass scale 
(ISesana et alJl2005b IVolonteri et al.ll2003l) . Following MMW, 
all masses, positions and softening lengths were rescaled by a 
factor H(z = 3)/ Hq, i.e. the ratio between the Hubble constant 
at z = 3 and its present-day value for a ACDM "concordance" 
cosmology (H = 70 km s" 1 Mpc" 1 , D, m = 0. 3, Q A = 0.7) 
The ha lo concentration was chosen according to Bulloc ket al.l 
(1200 ll) . c = 3. Satellite galaxies were initialized with the same 
structure, with the mass in each component being scaled down 
by q. The resulting SMBH pairs fall in the typical range 
of masses whose co alescences will be detectable with LISA 
(ISesana et al.l 120051) . We choose orbital parameters for the 
mergers that are common for merging halos in cosmological 
simulations (Benson 2005): the baricenters of the two galax- 
ies were placed at a distance equal to the sum of their virial 
radii and set on parabolic orbits with pericentric distances of 
20% the virial radius of the most massive halo. All mergers 
we considered were coplanar and prograde. 

All simulations wer e performed with GASOLINE, a 
TreeSPH N-body code dWadslev et all 120041) . We ran colli- 
sionless ("dry", with f g = 0) and gasdynamical ("wef) merg- 
ers with the same gas fraction in the primary and secondary 
galaxies, either / g = 0.1 or 0.3. In wet runs, atomic gas cool- 
in g was allowed ; star f ormation (SF) was treated according 
to IStinson et al.l (12006b . Gas particles are eligible to form 
stars if their density exceeds 0. 1 cm" 3 and their temperature 
drops below 1 .5 x 10 4 K, and the energy deposited by a Type- 
II supernova on the surrounding gas is 4 x 10 50 erg. With 
this choice of parameters our blast-wave feedback model was 
shown to produce realistic galaxies in cosmological simula- 
tions dGovernato et al.ll2007 ). A summary of our set of simu- 
lations is presented in Table Q] 

In each galaxy (except for a very high-resolution test, see 
13. p , we employed 10 6 particles for the halo, and, initially, 
2 x 10 5 star particles and 10 5 gas particles, when included. 
The force softening was 100 pc in our reference model, scaled 



TABLE 1 
Summary of Simulations 



q 


SF 


ft 


z 


BH final distance" 


0.25 


no 








2 - 4 kpc 


0.25 


yes 


0.1 





200 pc 


0.1 


no 





3 


1-6 kpc 


0.1 (hi-res) 


no 





3 


1 - 5 kpc 


0.1 


yes 


0.1 


3 


400 pc 


0.1 


yes 


0.3 


3 


70 pc 



11 When possible, estimates of pericenter and apoc- 
enter of the orbit of the lighter SMBHs inside the 
merger remnant are given. 

down by q 1 ^ in the satellites, and by H(z)/Hq in high-z runs, 
yielding a force resolution of ~ 20 pc in the satellite galaxy for 
q = 0.1 at z = 3. With such a high particle number, the masses 
of star particles in the satellite is an order of magnitude lower 
than Mbh, ensuring that SMBH dynamics is not affected by 
spurious two-body collisions. In what follows, we define two 
SMBHs as a "pair" if their relative orbit shrinks down to a 
separation equal to twice the softening. From these distan ces, 
a SMBH binary may form in ~ 1 Myr dMaver et al.ll2007l) . 

3. RESULTS 
3.1. Collisionless Mergers 

In collisionless runs, the satellite is not able to dissipate 
energy gained through tidal shocks at pericentric passages 
(iGnedin et al.lll999t iTaffoni et alJl2003h . For q = 0.25, dy- 
namical friction on the dark matter halo of the more massive 
galaxy is efficient, and the satellite sinks down to a few ^10 
kpc from the center after 3 orbits. At that point the central 
density of th e satellite has decreased considerably because of 
tidal heating (Kazantzidis et al. 2004). Its innermost region is 
then tidally disrupted, leaving the small SMBH at a distance 
of a few kiloparsecs. The dynamical friction timescale has 
now greatly increased, because the mass of the small "naked" 
SMBH is orders of magnitude lower than that of the satel- 
lite's core that once surrounded it. No pair is formed, and the 
smaller SMBH is left wandering a few kiloparsecs away from 
the center of the remnant (Fig. [T). We note that estimating 
correctly the dynamical friction timescale of the SMBH from 
the simulation is not trivial in this regime, because the dark 
matter component is still dynamically important at kiloparsec 
distances from the center of the remnant. Even at high res- 
olution, the mass of the dark matter particles of the primary 
galaxy is comparable with that of the SMBH, hence dynam- 
ical friction could be altered by discreteness effects. How- 
ever, the dynamical friction timescale needed for the "naked" 
SMBH to reach the center of t he remnant can be estimated us- 
ing Chandrasekhar's formula dColpi et alJI 19991) . and it turns 
out to be longer than a Hubble time. Hence we conclude that 
the two SMBHs will not form a pair. 

In the q = 0.1 case, dynamical friction is rather ineffective. 
The sinking time for the satellite is longer than a Hubble time 
for a z = merger, owing to the low initial mass of the satel - 
lite and to mass loss due to tidal stripping (Colpi et al.ll 19991) . 
However, since mergers are much more common at higher 
redshift, when orbital times are shorter by a factor H(z) /Ho, 
we performed a q = 0. 1 merger starting at z = 3 (see an 
epoch at which these SMBH pairs are predicted to be most 
typical. This is completed in ~ 2.5 Gyr. Similarly to the 
q = 0.25 case, a wandering SMBH is left at several kiloparsecs 
from the center of the primary (Fig. [T). In order to check that 
the tidal disruption of the core was not affected by numerical 
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FIG. 1 . — Separation of SMBHs as a function of time in four of our simu- 
lations. Upper panel: SMBH distance in q = 0.25 mergers; the thin and thick 
lines refer to the dry and wet cases, respectively. The inset shows the color- 
coded density of stars (left) and gas (right) for the wet case at t = 5.75 Gyr 
(marked in red on the curve); each image is 12 kpc on a side, and colors code 
the range 10~ 2 - 1 Mq pc~ 3 for stars, and 10~ 3 - 10 - ' Mq pc~ 3 for the gas. 
Lower panel: SMBH distance for q = 0.1, z = 3; the thin and thick line refer 
to the dry and /„ = 0.3 cases, respectively. The inset shows density maps at 
t = 1.35 Gyr for the f g = 0.3 case; images are 4 kpc on a side (color coding 
as in upper panel). 

heating, we ran the same merger with a 5 times higher mass 
resolution in the stellar component of both galaxies, and cor- 
respondingly higher force resolution; no significant difference 
in the SMBH orbit was found (Tab. [TJ. Therefore, even in this 
case the two SMBHs do not form a pair. 

3.2. Dissipational Mergers with Star Formation 

The presence of a star-forming gaseous component cru- 
cially affects the orbital decay of the SMBH via its dynamical 
response to tidal forces and torques. 

The orbits of dry and wet, q = 0.25 mergers differ only after 
the first couple of orbits (~ 6 Gyr, see Fig. [TJ. At second 
pericenter, tidal forces excite a strong bar instability in the 
satellite. Dissipation in the gas and torques exerted by the 
stellar bar onto the gas drive a gaseous inflow toward the cen- 
ter of the satellite (see inset in Fig. [TJ, increasing the central 
star formation rate by a factor of 3. Thus, the potential well of 
the satellite deepens, ensuring resilience of its central part to 
tidal stripping and shocks even when it plunges near the center 
of the primary. As a consequence, the small SMBH contin- 
ues to sink fast, because it remains embedded in the massive 
core of the satellite. A pair of SM BHs is formed in this case, 
confirming previous results (Kaza ntzidis & al.ll2005l) . 

In q = 0.1, z = 3 wet mergers, both / g = 0.1 and 0.3 were 
employed; the latter should be a more realistic assumption, 
since disk galaxies at z = 3 are believe d to have a higher gas 
mass fraction (e.g., [Franx et al.ll2008l) . In these cases, star 
formation and supernovae feedback affect the structure of the 
inte rstellar medium (ISM) in the disks quite dramatically (see 
also iGovernato et aT]|2007t) . The disks develop a clumpy and 
irregular multi-phase structure, and turbulent velocities of the 
gas become a significant fraction (30%) of the circular veloc- 
ity in this mass range (V v i r = 64 km s" 1 ). Star formation in 
the center of the satellite is enhanced compared to the same 
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FIG. 2. — Cumulative bound mass profiles of z = 3, q = 0. 1 satellites. Thin 
red lines refer to the fg = 0.3 satellite in isolation, thick blue lines to the 
fg = 0.3 merging satellite, and dashed green lines to the merging f g = 0.1. All 
data refer to the third apocenter, or equivalent time in isolation (t = 1.8 Gyr). 
Left panel: total bound mass profiles. The more gas-rich satellite devel- 
ops a higher concentration during the merger, compared to the other cases. 
Tidal truncation and gas removal cause a factor of ~ 2 difference in mass at 
2 kpc between isolation and merging cases. Right panel: bound mass in stars 
formed after the start of the simulation. The total amount of stars formed 
depends roughly only on the initial fg, but SF is more localized in the center 
during the merger because of tidal forces. The more gas-rich merging satellite 
undergoes the strongest central SF burst, developing a higher central density. 



model evolved in isolation, even though the star formation 
rate integrated over the whole galaxy is unchanged (see Fig. 
[2j The concentrated star formation and the turbulent motions 
of the gas stabilize the disks against bar instability (see inset 
in Fig. [TJ. In the absence of bar-driven torques, the strong 
gaseous inflow seen for higher-mass objects (q = 0.25) does 
not happen for q = 0.1. Yet, tidal torques drive some gas to- 
wards the center, explaining the enhanced central star forma- 
tion (Fig. [2j. During the first three orbits, the f g = 0.1 case 
preserves its initial central density owing to the mild mass 
inflow, rather than lowering it as in the collisionless case, 
while the f g = 0.3 satellite develops a steeper stellar cusp. 
Once the satellites go through the second pericentric passage, 
their ISM is prone to ram pressure stripping by the gas disk 
of th e primary galaxy, out side the ram pressure stripping ra- 
dius (Ma rcohni et"alfl2003h . Nearly 90% of the gas is swept 
away when they first enter the disk of the primary, while what 
remains is stripped during the next orbit: at t = 2 Gyr, the 
satellites have lost all their gas content, even in their cen- 
tral region. From this point onward, the satellite with initial 
f g = 0.3 is a cuspy, gas-poor object, subject to dynamical fric- 
tion in the stellar and gaseous background. Its sinking is rel- 
atively fast because the steeper stellar density profile implies 
a larger bound mass, enhancing dynamical friction relative to 
the dry merger c ase. Moreover, its r esponse to tidal shocks is 
nearly adiabatic dGnedin et al.ll 19991) . preserving it from tidal 
disruption. On the contrary, the satellite of the f g = 0. 1 run un- 
dergoes a slower decay because of the lower bound mass. It 
then experiences a higher number of tidal shocks at pericentric 
passages which further decrease its density, until its complete 
disruption. As a consequence, the f g = 0.3 merger leaves the 
lighter SMBH at 70 pc from the more massive one (Fig. [JJ 
Tab. [TJ in a gas-rich environment, where the dynamical fric- 
tion timescale for the SMBH to sink to the center, based on 
Chandrasekhar's formula, is very short (< 1 Gyr). Instead, in 
the / g = 0. 1 case the final distance of the SMBHs is ~ 400 pc; 
at such separations the dynamical friction timescale is of a 
few billion years. Hence the pairing will occur within one 
Hubble time in both cases, but will be considerably faster in 
the / g = 0.3 case. 
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4. DISCUSSION AND CONCLUSIONS 

Our results show that the formation of a SMBH pair in 
unequal-mass mergers of disk galaxies is very sensitive to 
the details of the physical processes involved. None of the 
collisionless cases we studied led to SMBH pairing: tidal 
shocks progressively lower the density in the satellite until 
it dissolves, leaving a wandering SMBH in the remnant. The 
inclusion of gas dynamics and SF changes significantly the 
outcome of the merger. For higher mass ratios (q = 0.25) at 
z = 0, bar instabilities funnel gas to the center of the satellite, 
steepening its potential well and allowing its survival to tidal 
disruption down to the center of the primary. Therefore, in 
this case the presence of a dissipative component is necessary 
and sufficient to pair the SMBHs at ~ 200 pc scales, creat- 
ing favorable conditions for the formation of a binary. The 
smaller satellites here considered (q = 0.1, z = 3) are more 
strongly affected by both internal SF and the gasdynamical 
interaction between their ISM and that of the primary galaxy. 
Torques in the early stages of the merger are funnelling the 
gas to the center less efficiently, due to the absence of a stel- 
lar bar and the stabilizing effect of turbulence. As a result, 
ram pressure strips away all of the ISM of the satellite. If 
satellites develop a higher central stellar density by rapidly 
converting their gas into stars before ram pressure removes 
it, they can retain enough bound mass to ensure the pairing 
of the two SMBHs. Gas-rich satellites (f s = 0.3) undergo a 
stronger burst of SF during the first orbits, and therefore meet 
this requirement better than f g = 0.1 satellites. Yet in both 
cases the central density of the cusp remains high enough to 
permit its survival, allowing the pairing of the two SMBHs 
within a Hubble time. 

The pairing of the two SMBHs takes less than 1 Gyr in the 
gas rich systems that should be common at z = 3. Therefore, 
if the MBH-A^buige relation approximately holds at z = 3 as in 
the local Universe, the galaxies here considered should lead 
to the fo rmation of representa tive SMBH pairs at such cosmic 
epochs (IVolonteri et a l. 2003). These pairs are also expected 
to contribute significa ntly to the high-z gr avitational wave sig- 
nal in the LISA band (ISesana et al. 2005). Since we show that 
gas-dynamical processes allow such an efficient pairing of the 
SMBHs, the results of this Letter strengthens the case for the 
observability of these coalescence events. On the other hand, 
we show that the timing between galaxy mergers and merg- 
ers of their SMBHs is sensitive to the gaseous content of the 
merging galaxies. Hence, SMBH coalescence events do not 
necessarily trace galaxy mergers directly. This will have im- 
portant implications on the interpretation of the LISA data 
stream. 

We note that the orbital evolution of the SMBH pairs, in 
the dynamical range considered here, has only a weak depen- 
dence on the masses of the two SMBHs. As shown in Figure 
12 the stellar mass enclosed inside two softening lengths from 
the center of our galaxy models (hence close to our resolution 
limit) already exceeds Mbh by more than an order of magni- 
tude. This is the effective mass that determines how quickly 
the SMBHs will sink. Therefore, lowering Mbh or increasing 
it by up to an order of magnitude would have no effect on sink- 
ing timescales before the disruption of the satellite. Instead, 
after disruption, the analytic estimate for the dynamical fric- 
tion timescale of the naked SMBH would change linearly with 
Mbh- Similarly, if nuclear star clusters with m asses ~ 10 Mbh 
( Wehner & Harrijl |2006t iFerrarese & aT]|2006l) were present 
around the SMBHs, their sinking timescales would still be 



longer than a Hubble time in our dry mergers, where the fi- 
nal SMBH separation exceeds 1 kpc, while the pairing would 
now occur in well below a Gyr in all our wet mergers. Hence, 
either a larger Mbh or the presence of a nuclear star cluster 
would enhance even further the difference between dry and 
wet mergers. 

Lastly, a general limitation of our simulations is that they 
lack gas accretion onto the SMBHs and associated energy 
feedback. Additional heating from the active SMBH should 
reduce the binding energy of the gas, making it more suscep- 
tible to stripping processes, and perhaps inhibiting the forma- 
tion of a steep stellar cusp. This would reduce the efficiency 
of the pairing process, but the effect will strongly depend on 
when the SMBH becomes active during the merger. Although 
it is unlikely that these effects will change the overall picture 
presented in this Letter, they will have to be explored in a 
forthcoming paper. 
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